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Abstract 

Reconnection of the self-generated magnetic fields in laser-plasma interaction was first investi- 
gated experimentally by Nilson et al. [Phys. Rev. Lett. 97, 255001 (2006)] by shining two laser 
pulses a distance apart on a solid target layer. An elongated current sheet (CS) was observed in 
the plasma between the two laser spots. In order to more closely model magnetotail reconnection, 
here two side-by-side thin target layers, instead of a single one, are used. It is found that at one 
end of the elongated CS a fan-like electron outflow region including three well-collimated electron 
jets appears. The (> 1 MeV) tail of the jet energy distribution exhibits a power-law scaling. The 
enhanced electron acceleration is attributed to the intense inductive electric field in the narrow 
electron dominated reconnection region, as well as additional acceleration as they are trapped in- 
side the rapidly moving plasmoid formed in and ejected from the CS. The ejection also induces a 
secondary CS. 

PACS numbers: 52.38.Fz, 52.35.Vd, 52.30.-q, 52.50.Jm 



Magnetic reconnection (MR) in plasmas is associated with explosive conversion of mag- 
netic energy into plasma kinetic and thermal energies. In the process, plasma is accelerated 
in and ejected from a thin region where reconnection takes place |l|-|5J]. The reconnection rate 



obtained from the experiments and solar observations is often 
predicted by the standard Sweet-Parker and related models 



bund to be larger than that 
-6| and has been attributed 



;o various collisionless and higher-dimensional effects, including Hall current and turbulence 



7H12J. It has also been shown that the reconnection rate is enhanced by the formation 
of secondary magnetic islands and plasmoid ejection from the region, which can be highly 
unstable if the Lundquist number S > 10 4 [13J|. These theoretical predictions are consistent 
with in situ observation of a secondary magnetic island near the center of the ion diffusion 



region in the near Earth magnetotail 14 1. It is well known that in laser-matter interaction 



meg agauss magnetic fields can be generated through the barotropic (Vn e x VT e ) mechanism 



15l Il6( |. Nilson et al. [17j were the first to use two approaching laser-created plasmas to 



model the reconnection geometry. Definitive proof of MR taking place is provided by the 
measurements of Nilson et al. |l7t] and Li et al. [18( on magnetic topology changes using 
time-resolved proton deflection technique, and by the observations of Nilson et al. [13] on the 
highly collimated bi-directional plasma jets 40° off the expected MR plane. In this Letter, 
we investigate the structure of collisionless reconnection of the self-generated magnetic fields 
generated in laser-target interaction produced plasmas by using two coplanar plane targets 
with a gap in between, so that the colliding magnetic fields and plasmas are not electrically 
pre-connected. Three well-collimated energetic (MeV) electron jets, or elongated electron 
diffusion regions (EDRs), inside and along the edges of a fan-like plasma outflow from the 
thin current sheet (CS) are observed. The energy gain is attributed to electron acceleration 
by the intense local inductive electric fields as well as inside a fast propagating plasmoid 
ejected from the elongated CS during the MR. The plasmoid also modifies the local magnetic 
field as it moves forward, producing a secondary CS. Our results agree quite well with that 
from numerical simulations and in situ observations. 

The experiments were performed at the SG-II laser facility. The general setup is similar 
to that in our earlier work [&]. Here, two 0.7 mm x 0.3 mm x 50 /zm thin Al foils separated 
by 150-240 /im are juxtaposed in the same (vertical) plane. Two heater laser beams are 
simultaneously (within ±25 ps) focused by F = 3 optics onto each foil. The focal spots, 
400 ± 25/im apart, have diameters < 150/xm FWHM. The incident laser intensity at each 



spot is ~ O(10 15 ) Wcni -2 , and the total laser energy on each target is ~ 450 joules. The 
main diagnostics for the density of the laser-produced plasma is a modified Nomarski in- 
terferometer of magnification 3.5. A 527 nm 150 ps FWHM Gaussian probe laser beam is 
directed perpendicular to the target plane at suitable intervals with respect to the heater 
beams. A polarization imaging system is used to measure the self-emission profile of the 
plasma, and a combination of reflective and bandpass (A ~ 532 ± 1 nm) interference filters 
is used to limit the detectable emissions. Three pinhole cameras with magnification 10 are 
used to monitor the x-ray emission from the front, back, and side of the plasma. At the top 
front of the target plane, an adjustable-range RAP crystal spectrometer of ~60 /zm spatial 
resolution is used to characterize the plasma composition. 

Figures QTa) and (b) show the interferometric images of the laser-generated plasma at 
t = to and t = to + 1.0 ns, respectively, where t = to corresponds to the instant when the 
tail-side half-maximum of the heater pulse reaches the target front. By symmetry the MR 



region is centered in the gap between the laser spots [171. Il8|. Figure [Tfc) shows the profile 
of linearly polarized 532 nm emission from the plasma, with the darkest shade representing 
the brightest emission. Fig. QJd) shows the upper half of the MR region as obtained from 
a particle-in-cell (PIC) simulation using the experimental parameters (details given below). 
The main features in Fig. [Tfc) are sketched and identified in Fig. QJe), where the solid (red) 
lines show the brightest emission regions. All the panels in Fig. [1] and the left panel in Fig. [2] 
show the MR, or the x — z plane, where x is along the two laser spots. By following an edge, 
marked F in Fig. HJa), the local plasma advection velocity is determined to be ~ 2000 km/s 
[lOj. The CS region (from XI to X2) between the two laser spots is of length ~ 350/im. Here 
are two characteristic regions: Rl points to a fan-like (with full-opening angle 26 rcc < 40°) 
region of irregular plasma motion, and R2 points to a relatively calm and uniform plasma 
region. Fig. QJb) shows that after 1.0 ns the fan-like region has expanded and the jets 
lengthened. In the emission image in Fig. [Tfc) one can clearly see the two bright (solid) 
lines marking the edges of the magnetic-separatrix region (see also Fig. [He)). At the center 
of the fan-like region, a third bright line corresponding to a well-collimated energetic-electron 
jet appears. The three electron jets have similar widths (~ 55/xm) in the optical images. The 
variance of the lengthes of the three jets in the time-resolved interferometric images indicates 
that the center jet is much faster than that on the two edges of the outflow fan, which has 
a bulk velocity v c f > 600 km/s, and that the center part of the CS has lengthened. Similar 



jet-like electron outflows have been observed in the magnetosheath and the correspondin 



numerical simulations 19j, as well as around the separatrices of the reconnection layer 20]. 
The upper half of the MR region in Fig. [T] involves more complex behavior. There is also 
a fan-like outflow of plasma from the CS. However, away from the primary CS there is a 
secondary CS with local MR. The secondary CS is connected to a flare loop at one end and 
a plasmoid at the other end, which in Fig. QTe) is encircled by a thick black dashed line. 
Within the latter one can see two bright areas (marked by red solid lines) that correspond to 
the leading and trailing ends of the plasmoid. One can also see that the trailing end of the 
plasmoid is connected to the secondary current sheet. A plasmoid can also be seen in the 
simulation result shown in Fig. QJd). It is born inside the primary CS and grows rapidly. Its 
rapid ejection stretches the reconnected magnetic field lines, creating a secondary CS and 
MR. In Fig. QJc) of the experimental results, the plasmoid's upward-transport is manifested 
by its motion-blurred emission images. However, the generation process of the plasmoid 
needs further more detailed studies. These identifiable features in Fig. [l](c) are also marked 
in Fig. []Te) and shall be discussed in more detail later. Our results thus support that of 



Shepherd et al. [2jJ and Uzdensky et al. [221 ]. and should be useful for analyzing abrupt 



solar and astrophysical events 23|, |24[. The top-bottom asymmetry of the MR region is 
inevitable in the case with plasmoid ejection because the CS is highly unstable. Any small 
local disturbance can lead to rapid plasmoid ejection. That is, the ejection is stochastic, as 



proposed by Uzdensky et al. [22] 



More detailed characteristics of the primary CS are obtained from the x-ray images and 
the emission spectra. Fig. [2] shows the x-ray pinhole image taken from the back of the 
Al targets. Intense x-ray emission comes from a narrow region of hot and dense plasma 
between the laser spots, with less intense x-ray emission comes from a fan-like region below 
it. The former region, of width ~ 85/im, can be taken to be the ion diffusion region. A Ti 
wire of diameter 50/im placed ~ 250/mi above and parallel to the target plane is used to 
determine the extent of the primary CS normal to the MR plane. As shown in Fig. [2J plasma 
ejection appears on the lower side without the Ti wire, but that on the upper side is mostly 
blocked by the latter. The extent of the CS is determined to be about 50 /im, indicating 
a quasi two-dimensional (2D) nature of the MR region 19|] . However, at longer times the 
plasma system becomes highly irregular and difficult to follow with our diagnostics. This 
may be due to onset of transverse instabilities that lead to a strongly turbulent state 12]. 



The spatial distributions of the K a resonance lines (He-a lines) of the helium-like Ti ions 
and their satellites, as shown in Fig. [3]^a), are used to determine the cross-sectional size of 
the energetic electron outflow [25j|. Fig. E^b) shows the spatially integrated line spectra. 
Although the Ti wire is also in the transverse expansion paths of the two laser-produced Al 
plasmas, the He— a lines of Ti plasmas are only present in the midplane region (between 
80 — 150/im). Accordingly, together with the evidence given by the bright spot created by 
the central jet on the sampling sphere, as can be seen in the lower parts of Figs. Hfb), (c), and 
(e), one can conclude that a well-collimated energetic electron jet is formed. The uniform 
continuum emission on both sides of the Ti He-a lines is from the Al plasma. Simulations by 
Shay et al. suggest that the structure of the fan-like MR plasma outflow and well-collimated 
energetic electron jets can be attributed to competition and balance of the inductive electric 
field, the electron inertia, the Lorentz force, and divergence of the momentum flux [19] . 
However, more comprehensive laboratory and in situ magnetotail measurements are still 
needed in order to more definitively understand the detailed physics of the process |19|, l20| . 

Figure EIc) shows simultaneously the measured He-a and Ly-a line emissions of the Al 
plasma, and Figs. EJ^d) shows the space-resolved line spectra. The x-ray spectra of the 
Al plasma at the free-expansion edge, namely along P3 in Fig. Etc), shows a clear He-a 
line and a weak Ly-a line. From non-LTE calculations, we found that this plasma region 
has temperature ~ 440 eV and density ~ 2.5 x 10 19 cm~ 3 , which are consistent with the 
interferometric measurements. The x-ray spectra from the plasma at the center, namely 
along PI in Fig. EJ^c), clearly shows a saturated He-a line and a strong Ly-a line, as well as 
several satellite lines from Li-like ions, indicating an electron temperature of ~ 570 eV and 
density of ~ 4.8 x 10 19 cm~ 3 , the latter about double that at the free-expansion edge (i.e., 
along P3). Fig. Ele) shows the electron temperature and density from the experiments and 
hydrodynamic calculations. However, the Ti plasma produced by impact of the energetic jet 
electrons with the Ti wire has a temperature of ~2300±100 eV and density of 3.2±0.5xl0 23 
cm -3 , which are both considerably higher than that in the diffusion region. This is due to 
the fact that the collimated jet electrons have speeds much higher than the speed ~ 600 
km/s obtained from the intense linearly polarized 532 nm synchrotron radiation imaging of 
the plasma profile discussed above. 

The electron and ion inertial lengths as calculated are d e = c/u pe = 0.77— 1.06/im and di = 
c/ojpi ~ 54 — 75/im, respectively, where c is the light speed, and cu pe (uj p i) is the electron(ion) 



plasma frequency. The electron density is n e ~ 2.5 — 4.8 x 10 19 cm" 3 , and the average ionic 
charge is Z ~ 10. Accordingly, the width ~ 85/zm of the defined ion diffusion region in 
Fig. H] is ~ 1 — 1.5 times the ion inertial length. The widths of the EDRs are obtained by 
two methods, which are sensitive to the photon energy ranges. The pinhole x-ray image 
at >1.5keV gives the widths of the EDRs as 7 and 15 /mi as in Fig.2, corresponding to 
~ 10 — 15d e , respectively. On the other hand, the polarized emission image at 532 nm 
gives the widths of the three EDRs as ~ 55/mi. The difference in the EDR widths from 
the two measurement methods indicates that, as expected, the center regions of the EDRs 
are of higher temperature. The electron mean-free-path is calculated as \ em f p ~ 114/im, 
which is larger than the measured current layer width. For the experimental parameters, we 
have the ratio between the ion skin depth and the Sweet- Parker width djld sp ~ 10 > 1, so 
that the Hall effect should be the main effect for the fast reconnection [26(, and the binary 
collision is thus negligible in the CS, as also assumed for our PIC simulations. With detailed 
characteristics of the MR plasma region, and the separatrix magnetic field of So ~ 3.75 MG 
inferred later, we now can the give the upper limit of the Alfven speed as Va ~ 720 km/s. 
The measured velocity ~600km/s of the two side EDRs in Fig. [1] is near the calculated value 
of Va- However, the central EDR that starts later than the two side ones has clearly a higher 
propagation speed (since it can almost catch up with the latter), perhaps super-Alfvenic as 
predicted by PIC simulations. 

We have carried out 2D PIC simulations of the main MR region using parameters close to 
that of our experiments. However, a large (mj/m e =100) electron mass and a relatively low 
(c/Va = 30) light speed are used [27| . The initial state is a Harris current sheet of width di 
in a L x x L z = 12.8di x 102.4g?j simulation box. A guide field B y ~ 0.2S , where B z is the 
strength of the sheared magnetic field, is added to model the experimental magnetic topology. 
Initially, the plasma temperature is T i0 = T e0 = 0.011m e c 2 , and the plasma beta is f3 = 0.2. 
Fig. [T](d) shows a simulated current density distribution j y (color coded) in the upper halfoi 
the MR region, together with several typical field lines (black lines). Here one also finds a 
plasmoid, together with the deformed magnetic field lines and a secondary CS, as observed in 
the experiments. However, it should be noted that our simulation does not include the initial 
laser-plasma interaction stage since that would be beyond our computational resources. In 
fact, it only models the thin central MR region, which should be qualitatively similar for 
most MRs [7| , and the initial values of the plasma and field parameters have also been taken 



from the already evolving system in the experiment. 

The energetic jet electrons from the diffusion region are also diagnosed by using a mag- 

The jet electron energy can reach MeVs and obeys a power- 



netic spectrometer 



law (of index —6.5) scaling at the high-energy tail for electrons up to ~ 2 MeV (the 
limit of our measuring equipment). We have also estimated the energy of the relativis- 
ts electrons from the synchrotron radiation (see Fig. C](c)). For relativistic electrons gy- 
rating in the MR and Hall magnetic fields, the maximum electron energy £ max is related 
to the maximum measured wavelength A max of the linearly polarized synchrotron radia- 
tion by £ max = m e c 2 (2.5 x 10 _2 /A max -B ) 1 ^ 2 , and the corresponding electron gyroradius is 
r g = 1.70 X 10 3 A/(£max/0.511) 2 — 1.0/£>o- Assuming that the width of the electron diffusion 
region is of the order of r g ~ 55/zm,and the maximum wavelength is A max = 532nm, we 
obtain 5.2 MeV as the measured energy of the relativistic electrons and B ~ 375 Tesla as 
the strength of the reconnected magnetic field. Existing numerical simulation and analytical 



estimate indicate that electrons can be accelerated to ultrahigh energies by the e 
produced in MR, which explains the good collimation of energetic electrons [27|, 



ectric fields 



2J-l3l|. For 



electrons to achieve MeV energy in ~ 100 ps (the tolerant time for Fig. Q](c)) to prevent 



from being blurred by plasma expansion; also see 17|, ll8|), the electric field strength should 
be ~ 10 8 V/m. According to Faraday's law, such a reconnection electric field corresponds 
to the annihilation of several tens Tesla of magnetic fields in ~ 0.1 ns, or a dimensionless re- 
connection rate [19( 0.1 - 0.3, which is consistent with the Hall MR rate ~ tan(# rec ) ~ 0.25 
- 0.35 obtained from measuring the local magnetic exhaust wedge angle 28 rec in Fig. [TJ 
The corresponding acceleration length should then be ~ 1cm. That is, the pre-accelerated 
electrons trapped by the rapidly moving plasmoid have bounced many times before they 



gain sufficient energy to escape 32| ) . Accordingly, the MR electric field and plasmoid trap- 
ping, accounting for the good collimation and the power-law scaling, respectively, are both 
responsible for the observed characteristics of the high-energy electrons. 

In summary, two suitably juxtaposed Al foil targets are used to produce in the laser 
generated plasma the configuration of MR in the presence of a guide field. It is found that 
in the MR a plasmoid is generated and ejected from the CS. As it rapidly propagates away 
it deforms the reconnected magnetic field and generates a secondary CS as well as flare 
loops. The primary MR also involves fan-like plasma outflows, in which appears intense 
well-collimated electron jets containing MeV electrons. The results are consistent with the 

8 



theoretical prediction of anomalous plasmoid generation by unstable CS, often invoked to 
interpret solar flares. Simulations show that the guide field can indeed lead to the fast 
formation of the plasmoid and the secondary CS. The results here also demonstrate that 
suitably designed laser-plasma interactions can be effective and economical for laboratory 
modeling high-energy-density astronomical and astrophysical phenomena. 
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FIG. 1: (Color) Interferometric image at t = to ns (a) and t = to + 1 ns (b), where to is set at the 
half maximum of the laser-pulse tail. The electron diffusion region is between XI and X2. F marks 
the edge of the advecting plasma. PI and P2 (1.16 mm and 1.75 mm from X2, respectively) mark 
the location of a front of the expanding fan-like region at the two times, from which an outflow 
speed of > 600 km/s is estimated. The angle of the fan- like plasma region is 35° - 40°. (c) Linearly 
polarized self-emission image of the plasma, (d) PIC simulation result showing the current density, 
typical field lines, the secondary CS, and a plasmoid. (e) Schematic illustration of (c). The short 
thick (red) arrows indicate the plasma inflow directions for the primary and secondary CS. 
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FIG. 2: A typical x-ray pinhole image, with superimposed schematics of the experimental setup, 
and the corresponding intensity curves at two locations. The x-ray emission is strongest in the 
electron diffusion region at the center, and is of width ~ 85/iin. The two central emission peaks in 
the L2 profile correspond to the two edge electron jets at the start of the fan-like region. 



Laser foci on 



Lwer 1 



e50(.un 
_ Ti wire 

-F2 t2±. hot jets 11% „ 

U Imt 5p 




47CO 4725 4750 477S 157& 1B9C 1EC5 172S 
Pnoto energy feVj Photo enef^j (eV) 



■200 -100 100 200 
Position pfli} 



FIG. 3: (Color) Spatially resolved emission of K a lines from helium-like ions and the satellite lines 
from (a) the Ti wire, and (c) the Al targets, (b) and (d) Typical line spectra from Ti and Al, 
respectively, together with the results calculated using a non-LTE model, (e) Local Al plasma 
electron density and temperature estimated from (d). The solid and dashed lines show the plasma 
electron density and temperature, respectively, obtained from a hydrodynamic code. The sketch 
shows the relevant experimental components. 
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